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Coupled  CFD/FE  heat  transfer  /  heat  flux  calculations  for  material  temperature  prediction 
with  independent  grids  expose  the  need  to  exchange  boundary  condition  information  between 
the  fluid  and  the  solid  domain.  The  author  introduces  a  method  of  storing  the  boundary  con¬ 
dition  information  together  with  the  geometry  information,  where  boundary  conditions  are 
represented  by  a  set  of  geometric  entities.  This  method  of  “virtual  boundary  conditions”  sim¬ 
plifies  some  aspects  of  the  coupled  calculations  by  introducing  a  clear  grid  interface  and  en¬ 
ables  a  master  model  approach  from  the  geometry  generation  (CAD)  down  to  the  analytical 
disciplines  (CAE).  The  detachment  of  the  boundary  condition  information  from  the  computa¬ 
tional  grid  also  exposes  advantages  when  considering  different  degrees  of  geometric  com¬ 
plexity  (2D/3D). 


Introduction 

Modem  jet  engines  with  their  highly  loaded  components  expose  a  need  to  use  advanced  three-dimensional 
analytical  tools  to  predict  correct  life  times  during  design.  The  progress  in  computer  aided  design  (CAD)  and 
computer  aided  engineering  (CAE)  has  enabled  the  use  of  three  dimensional  models  and  computations  for 
almost  every  part  of  a  jet  engine.  Coupled  heat  transfer  /  heat  flux  calculations  represent  a  tool  for  very  accu¬ 
rate  prediction  of  temperatures  in  structural  elements.  They  calculate  the  amount  of  heat  taken  into  or  from 
the  structural  element  by  the  surrounding  gas  flows  by  modeling  the  gas  flows  themselves  and  the  heat  flux 
inside  of  the  structure.  Nevertheless,  these  calculations  require  extensive  exchange  of  information  between 
the  different  computational  domains,  which  makes  this  information  exchange  a  key  problem  to  overcome. 
Even  though  it  is  possible  today  to  automatically  create  grids  from  CAD  models,  it  is  still  very  uncommon  to 
also  represent  boundary  conditions  or  calculation  results  inside  the  CAD  model.  Therefore,  the  idea  of  a 
master  model  that  drives  the  entire  design  up  to  manufacturing  is  currently  excluding  the  analytical  disci¬ 
plines.  In  this  paper,  the  author  introduces  a  method  of  storing  and  exchanging  boundary  conditions  for  cal¬ 
culations  together  with  the  CAD  geometry  in  order  to  simplify  the  design  process  and  coupled  calculations 
in  order  to  come  a  step  closer  to  generating  a  complete  master  model. 

The  concept  of  “virtual”  boundary  condition  objects 

A  typical  design  process  starts  from  simple  ID  and  2D  geometries  and  evolves  through  2.5D  until  the  com¬ 
plete  3D  model  is  designed.  The  analytical  disciplines  accompany  this  process  by  increasingly  complex  cal¬ 
culations  where  the  boundary  conditions  stay  always  attached  to  the  used  grids  and  not  to  the  underlying 
geometric  model.  The  idea  of  attaching  the  boundary  conditions  to  the  geometry  rather  than  the  grid  makes 
use  of  the  fact  that  almost  every  software  describes  curves  and  surfaces  in  a  parametric  fashion  -  mostly  B- 
splines.  A  boundary  condition  can  now  be  regarded  as  a  geometric  element  itself,  which  is  called  a  “virtual” 
boundary  condition  that  is  associated  with  an  element  of  the  geometric  model.  Picture  1  depicts  how  a  tem¬ 
perature  boundary  condition  T  for  a  model  curve  is  stored  as  a  “virtual”  curve  with  identical  parameteriza¬ 
tion  in  the  plane  y=0.  For  easy  access,  the  x-coordinate  of  the  virtual  curve  is  chosen  to  be  identical  with  the 
u-parameter  and  the  z-coordinate  corresponds  to  the  value  of  the  boundary  condition. 
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Geometry 

Spline  curve  in  3D  xyz-space 
with  curve  parameter  u 


Pi 


Virtual  BC 

Virtual  boundary  condition  in  the 
y=0  plane.  The  boundary  condition  T 
is  mapped  as  a  function  of  u  for  each 
point  of  the  curve. 


Coordinate: 
Curve  parameter: 


u=0  u=u1  u=umax 


Figure  1  .  Schematic  representation  of  a  geometric  entity  and  a  corresponding  virtual  boundary  condition  element  in 
2D  space. 


For  each  possible  point  pi  on  the  geometric  curve,  a  boundary  condition  value  can  now  be  assigned  by  re¬ 
questing  the  curve  parameter  iq  of  that  point  and  reading  the  z-coordinate  of  the  virtual  curve  at  x=Ui.  With 
the  storage  of  boundary  condition  information  in  this  manner,  a  new  curve  is  necessary  for  each  boundary 
condition  (e.g.  temperature,  pressure  etc.)  that  should  be  associated  with  the  original  curve.  This  leads  to  a 
structure,  in  which  each  curve  uses  a  set  of  corresponding  virtual  boundary  condition  curves  to  represent  the 
entire  boundary  condition  information.  The  same  schematic  can  also  be  applied  to  surfaces  that  have  a  (u,v) 
parameterization.  For  each  surface  in  the  model,  a  set  of  virtual  surfaces  can  be  created  to  store  the  boundary 
condition  information.  The  u-parameter  of  these  virtual  surfaces  corresponds  to  the  x-coordinate  and  the  v- 
parameter  to  the  y-coordinate,  while  the  z-coordinate  holds  the  value  of  the  boundary  condition  respectively, 
as  can  be  seen  in  figure  2. 


Geometry 


Virtual  BC 

Spline  surface  from  a  CAD  model 
w\\h(u,v)  parameterization 


Virtual  boundary  condition  surface  with 
corresponding  (u;v)  parameterization  for 


Figure  2:  Representation  of  a  surface  and  a  corresponding  virtual  boundary  condition  element  in  three 
dimensional  space. 
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u-parameter 


Figure  3:  Creation  of  2.5  dimensional  geometric  objects  from  2D  curx’es  via  lofting.  The  virtual  boundary  condition 
objects  are  lofted  as  well  to  create  a  corresponding  2.5D  boundary  condition. 


This  way  of  storing  the  boundary  condition  information  has  advantages  in  the  progress  of  a  typical  design 
process.  Since  2  dimensional  models  and  analysis  characterize  the  early  stages  of  a  new  design,  this  infor¬ 
mation  can  later  be  used  to  switch  to  higher  dimensional  models  as  can  be  seen  in  figure  3. 

In  this  example  of  an  airfoil  design,  the  shape  of  the  airfoil  is  represented  by  three  curves  C1-C3  with  corre¬ 
sponding  virtual  boundary  condition  curves  T1-T3.  The  generation  of  2.5  dimensional  geometry  is  earned  out 
by  lofting  a  surface  through  the  curves  C1-C3.  The  boundary  condition  for  this  new  surface  is  created  ac¬ 
cordingly,  by  lofting  the  virtual  boundary  condition  curves  TrT3  to  create  a  virtual  boundary  condition  sur¬ 
face.  The  boundary  condition  information  can  then  be  used  as  a  stalling  point  for  the  2.5D  analysis. 

Accessing  the  information  of  virtual  boundary  conditions 

Once  the  boundary  condition  information  is  stored  in  the  “virtual”  form,  that  means  as  geometric  entities,  it 
can  easily  be  accessed  by  any  program  simply  by  dropping  the  (x,y,z)  coordinates  of  each  grid  node  onto  the 


Geometry  +  Virtual  Boundary  Condition 


Mesh  creation  Query  for  boundary  conditions  for  all  nodes  that  lie  on 

from  geometry  the  bounding  geometry 


Figure  4:  Extraction  of  boundary  condition  values  from  the  virtual  boundary  conditions.  A  point  on  a  mesh  that  was 
created  from  an  existing  geometry  (1.)  is  dropped  onto  the  geometry  (2.)  and  boundary  condition  values  are 
extracted  from  the  virtual  objects  (3.).  These  values  are  then  imposed  onto  the  mesh  (4.). 
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surface  and  obtaining  the  (u,v)  parameters  for  this  location.  With  these  parameters,  the  z-coordinate  value  of 
the  “virtual”  boundary  condition  surface  can  be  read  out  and  returned  to  the  analytical  program  as  boundary 
condition  value  at  this  node  as  visualized  in  figure  4.  The  close  interaction  between  CAD,  where  the  geome¬ 
try  and  the  boundary  conditions  are  stored,  as  well  as  CAE  becomes  obvious  with  this  approach. 

The  computational  power  required  for  the  determination  of  the  u  and  v  parameters  as  well  as  the  extraction 
of  the  boundary  condition  at  any  given  point  is  small  because  highly  optimized  subroutines  exist  for  these 
tasks  in  virtually  every  standard  software  package  that  uses  spline  surfaces. 


Application  for  coupled  calculations 


Coupled  calculations  usually  incorporate  different  solvers 
for  the  domains  to  be  investigated.  A  common  application 
is  the  coupling  of  CFD  solvers  for  fluids  and  FE  solvers 
for  solid  structures.  Each  solver  uses  a  grid  to  resolve  its 
computational  domain  and  the  exchange  of  information 
has  to  be  managed  by  interpolation  between  grids.  Figure 
5  depicts  that  situation,  showing  the  different  grids  and 
their  common  border  curve  where  the  grids  need  to  ex¬ 
change  information.  The  use  of  virtual  boundary  condi¬ 
tion  objects  does  not  eliminate  the  need  for  interpolation, 
but  it  significantly  reduces  the  complexity  of  the  problem 
by  creating  a  clear  interface  for  information  exchange. 
Nevertheless,  for  coupled  calculations  the  situation  is 
more  complex  since  the  virtual  boundary  condition  ob¬ 
jects  have  to  be  recreated  from  the  results  of  each  single 
computation.  The  results  are  typically  attached  to  the 
computational  grid  and  have  no  connection  to  the  under¬ 
lying  geometry.  Therefore,  this  information  has  to  be 
stored  in  an  independent  file  from  where  it  can  be  ac¬ 
cessed  to  recreate  a  boundary  condition  object  after  the 
calculation. 


Mesh  for 
calculation  1 


Common  border 
curve 


Figure  5:  Schematic  process  of  coupled  heat 
transfer/heat  flux  calculations  using  independent 
meshes  for  the  fluid  and  solid  domain.  The  virtual 
boundary  condition  elements  are  used  to  ex¬ 
change  information  at  the  common  border. 


The  general  process  used  for  coupled  calculations  is  explained  in  figure  6.  The  first  step  (1)  creates  the  grid 
for  one  computational  domain  (e.g.  the  solid  region)  from  the  geometry  where  the  (u,v)  parameter  informa¬ 


tion  for  each  node  is  stored  externally  (2).  This  information 
together  with  the  virtual  boundary  conditions  leads  to  the 
starting  grid  with  imposed  boundary  condition  values  at 
each  node  (3)  from  where  the  solver  calculates  a  first  solu¬ 
tion  (4)-(5).  These  results  are  then  re-transformed  into  vir¬ 
tual  boundary  condition  objects,  again  with  the  use  of  the 
(u;v)  parameter  data  (6).  Then,  the  whole  process  is  re¬ 
peated  for  the  other  computational  domain  (e.g.  the  fluid 
region)  until  a  converged  solution  is  obtained.  It  is  not  nec¬ 
essary,  that  all  boundary  conditions  are  rebuilt  after  each 
run,  sometimes  it  is  even  impossible  because  the  boundary 
condition  is  input  and  not  a  result.  For  the  CFD/FE  calcu¬ 
lation,  for  example,  the  boundary  conditions  for  the  fluid 
domain  can  be  constant  temperature  from  which  the  wall 
heat  flux  is  computed  and  added  to  the  list  of  boundary 
condition  objects.  The  solid  domain  calculation  uses  this 


boundary  condition  to  calculate  a  material  temperature. 


which  then  modifies  the  existing  temperature  boundary  Figure  6:  Complete  process  of  probing  for  the 


condition  for  the  fluid  domain.  By  repeating  this  scheme,  a  11  se  and  recreation  of  boundary  conditions 
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coupled  heat  flux  -  heat  transfer  calculation  can  be  performed.  The  obvious  difficulty  with  this  approach  is 
the  recreation  of  the  virtual  boundary  condition  object  from  the  discrete  (node  or  element  centered)  values  of 
the  calculation  results.  This  represents  a  best-fit  problem  for  curves  or  surfaces  from  a  cloud  of  points.  With 
block  structured  meshes,  the  best  fit  of  a  surface  is  still  fairly  simple  because  the  nodes  are  ordered.  For  the 
general  case  of  unstructured  meshes,  the  interpolation  of  a  surface  through  the  resulting  points  is  more  com¬ 
plex  as  can  be  seen  in  figure7. 


Interpolated  BC  Surface 


Figure  7:  Interpolation  of  a  new  virtual  boundary  condition  surface  from  result  point  data 


Surface  fitting  from  calculated  results 

There  are  many  good  algorithms  for  surface  fitting  through  irregular  data  points  and  a  good  overview  is 
given  in  [8],  For  the  sake  of  simplicity  in  programming  and  its  good  performance  with  large  numbers  of 
fairly  regularly  distributed  points,  a  cell  oriented  interpolation  method  similar  to  the  modified  Shepard’s  al¬ 
gorithm  developed  by  RENKA  [6]  with  a  regular  interpolation  grid  was  chosen.  The  algorithm  creates  a  rec¬ 
tangular  grid  of  interpolation  points  in  the  point  cloud  through  which  the  surface  can  then  be  fitted  easily  as 
show  in  figure  8. 
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Grid  of  Interpol.  Points]  [Interpolated  Surface] 
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Figure  8:  Creation  of  a  regular  grid  of  interpolation  points  used  to  create  a  surface  from  a  point  cloud. 


The  first  step  to  create  an  interpolation  point  is  to  find  its  closest  neighbors  and  then  use  the  distance- 
weighted  interpolation  with  these  points  to  calculate  the  needed  values  RiS 


RiS=- 


M  d 

m  =  l  “m 

M  , 

m=  1  Ct 


(SYB)  34-6 


with: 

R:  Value  of  boundary  condition 

i:  Index  of  boundary  condition  (in  the  case  of  multiple  boundary  conditions,  e.g.  pressure  and 
temperature) 

S:  Interpolation  point  with  parameters  (us,vs) 

m:  Index  denominating  all  neighboring  points  for  interpolation  (up  to  M) 

d:  Distance  between  point  Pm  and  interpolation  point  S 

With  the  cell  oriented  algorithm,  the  number  of  comparisons  needed  to  find  the  closest  neighbors  decreases 
dramatically,  because  only  points  in  concentrical  rings  of  neighboring  cells  have  to  be  considered  as  shown 
in  figure  9.  Depending  on  the  number  M 
of  data  points  considered  for  the  interpo¬ 
lation  of  each  grid  point  and  the  cell  size, 
a  numerical  smoothing  can  be  imposed  to 
the  resulting  surfaces,  leading  to  a  stabi¬ 
lization  of  the  coupled  computing  algo¬ 
rithm  described  above.  For  the  test  cal¬ 
culations  presented  in  this  paper,  eight 
surrounding  points  for  each  interpolation 
point  produced  good  results.  In  most 
cases,  the  amount  of  computing  power  is 
reduced  significantly  when  compared  to 
the  direct  interpolation  from  one  grid  to 
another,  because  only  relatively  few  in¬ 
terpolation  points  have  to  be  determined 
to  span  a  surface  over  the  resulting  data. 

A  coupled  heat  flux  /  heat  transfer  calculation  with  virtual  boundary  conditions 

As  described  above,  the  coupling  of  heat  flux  /  heat  transfer  calculations  is  used  to  demonstrate  the  method 
of  “virtual”  boundary  conditions  because  of  the  need  to  exchange  information  over  boundary  conditions 
between  the  fluid  and  solid  domain.  The  geometry  used  to  perform  the  calculations  was  experimentally  in¬ 
vestigated  by  BERG  [1],  Elfert  [3]  and  HEIN[10].  It  consists  of  a  tube  with  a  length  of  20  hydraulic  diame¬ 
ters  and  a  smoothly  rounded  inlet  area  as  shown  in  figure  10.  The  setup  was  designed  to  investigate  the  in¬ 
fluence  of  rotation  to  local  Sherwood/Nusselt  number  distributions.  It  was  later  used  to  calibrate  several 
CFD  calculations  ([5]  [2]  [7]).  In  extension  of  the  experiment,  a  computation  makes  it  possible  to  investigate 
the  influence  of  wall  thickness  and  external  heat  transfer. 
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Figure  9:  A  cell  oriented  method  reduces  the  amount  of 
comparisons  needed  to  determine  the  next  neighbors  by 
searching  in  cells  around  the  interpolation  point. 


Figure  10:  Schematic  layout  of  the  experimental  setup  used  by  Berg  [1  ]  and  Elfert  [3], 
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Therefore,  the  wall  is  modeled  as  solid  domain  for  the  heat  flux  calculations  and  the  inside  of  the  tube  is 
modeled  as  the  fluid  domain.  Figure  1 1  shows  a  coarse  version  of  the  grids  used  for  the  calculations  with  a 
modified  logarithmic  wall  function  for  heat  transfer  on  the  inside  of  the  pipe  and  constant  heat  transfer  from 
the  environment  of  the  tube. 

The  coupling  was  established  using  the  method  introduced  by  FlESELHAUS  [9]  in  order  to  guarantee  a  stable 
algorithm.  This  method  calculates  artificial  heat  transfer  coefficients  with  local  temperature  differences  at 
each  node  of  the  wall  and  uses  these  as  input  for  the  heat  flux  calculation.  After  the  finish  of  the  heat  flux 
calculation,  the  resulting  solid  material  temperature  is  used  again  as  boundary  condition  for  the  fluid  flow. 
The  entire  setup  consists  of  two  independent  grids  derived  from  several  surfaces,  each  equipped  with  a  set  of 
“virtual”  boundary  conditions  to  store  and  exchange  the  necessary  data  between  these  grids. 


Figure  11:  Computational  domains  used  for  coupled  calculation  (solid  domain  on  the  left  and  fluid  domain  on  the 
right) 

The  geometry  was  created  using  the  program  I  CAD  from  Knowledge  Technologies  International  (KTI) 
which  provides  a  geometry  kernel  in  combination  with  an  object  oriented  programming  language.  The  con¬ 
cept  of  “virtual”  boundary  conditions  was  employed  by  creation  of  a  new  surface  class  incorporating  a  list  of 
supplementary  surfaces  to  represent  the  boundary  conditions  as  attribute.  Methods  to  query  boundary  condi¬ 
tion  values  and  to  recreate  boundary  condition  surfaces  from  a  cloud  of  points  were  employed  and  linked  to 
that  surface  class.  The  CFD  computations  were  carried  out  using  CFX  TASCflow  V2 . 10  with  a  standard 
k-£  turbulence  model  as  well  as  modified  logarithmic  wall  functions  [4],  The  known  shortcomings  of  this 
approach  for  the  prediction  of  accurate  heat  transfer  were  accepted  because  the  focus  of  the  calculations  was 
the  proof  of  concept  for  the  boundary  condition  representation.  The  heat  flux  calculations  for  the  solid  do¬ 
main  were  carried  out  using  the  MSC  Qtran  software. 


Figure  12:  Velocity  vectors  showing  the  secondary  flow  field  at  the  end  of  the  pipe  (left)  and  comparison  with  the 
experiments  ofElfert  [3], 
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Figure  12  shows  the  visualization  of  the  secondary  flow  field  at  the  end  of  the  tube  at  x/dh=20  for  a  case  with 
Re=45000  and  Ro=0.044  as  well  as  a  comparison  of  the  main  flow  velocity  u  in  the  plane  of  rotation  (sym¬ 
metry)  and  the  secondary  flow  v  perpendicular  to  the  plane  of  rotation.  The  comparison  exposes  a  good 
match  of  the  flow  field  to  the  experiments  conducted  by  ELFERT  and  FlEIN. 
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Figure  13:  Material  temperatures  in  the  pipe  wall  for  different  wall  thickness 


With  the  coupling  algorithm  by  FlESELHAUS,  a  material  temperature  in  the  wall  of  the  pipe  was  calculated 
with  a  constant  environment  temperature  TE=800K  and  constant  heat  transfer  coefficient  h=400  on  the  out¬ 
side  of  the  pipe.  The  fluid  inlet  temperature  Tin  was  290K  and  the  fluid  heated  up  to  about  400K  at  the  out¬ 
let.  The  material  temperature  for  a  wall  thickness  of  dw/dh=  1 .5  and  dw/dh=  1.0625  is  shown  in  figure  13. 

With  this  information,  the  Nusselt  number  distribution  on  the  inside  of  the  pipe  can  be  calculated.  This  dis¬ 
tribution  shows  some  discrepancies  with  the  results  obtained  by  BERG  which  are  easily  explainable  by  the 
use  of  the  simple  k-£  turbulence  model  and  wall  functions.  Nevertheless,  the  general  trend  as  displayed  in 
figure  14  is  in  good  accordance  with  the  experiments  published  in  [1], 
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Figure  14:  Nusselt  number  distribution  on  the  inside  of  the  pipe  wall. 
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Conclusions 

The  method  of  virtual  boundary  condition  objects  is  implemented  and  works  well  for  coupled  heat  transfer  / 
heat  flux  computations.  With  the  geometric  representation  of  boundary  conditions,  a  close  connection  be¬ 
tween  CAD  and  CAE  is  possible  whereas  the  “virtual”  objects  act  as  interfaces.  The  boundary  condition  is 
stored  in  one  central  place  together  with  the  geometry  and  can  be  accessed  by  any  grid  generator  to  impose 
values  onto  the  nodes  generated.  Therefore,  the  boundary  condition  becomes  independent  of  the  used  grid 
type  which  enables  easy  comparison  of  results  obtained  with  different  meshes.  The  interpolation  time  can  be 
reduced  because  only  relatively  few  points  are  needed  to  describe  the  boundary  condition  surface.  The 


(SYB)  34-9 


boundary  condition  follows  the  geometry  generation  when  the  step  from  2D  to  3D  is  carried  out.  The  inter¬ 
polation  from  resulting  nodes  to  a  smooth  surface  works  well  and  can  be  carried  out  by  effective  interpola¬ 
tion  schemes.  Together  with  modern  CAD  systems  and  grid  generators,  the  representation  of  boundary  con¬ 
ditions  as  “virtual”  surface  objects  extends  the  master  model  philosophy  into  the  analytical  disciplines. 
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Paper  Number:  34 

Name  of  Discusser:  J.  J.  Me  Guirk,  Loughborough  University,  England 

Question: 

My  question  concerns  the  interpolation  procedure  for  generating  new  virtual  boundary  conditions.  Before 
the  interpolation,  the  temperature  filled  in  the  point  cloud  will  imply  a  certain  total  heat  flux  over  a  surface 
from  (say)  the  fluid  flow  into  the  solid  material. 

How  do  you  know  that  that  same  heat  flux  is  represented  by  the  interpolated  surface? 

If  not  then  the  heat  flow  is  not  concerned  between  fluid  and  solid  regions,  the  interpolation  procedure  may 
create  or  distroy  heat  flux? 

Answer: 

The  conservation  is  necessary  in  any  kind  of  interpolation  Scaling  of  the  virtual  boundary  condition  surface 
can  be  one  way  of  adjustment  to  account  for  this. 

A  comparison  of  the  calculated  total  heat  flux  and  the  generated  boundary  condition  surface  has  to  be  carried 
out  to  determine  the  scaling  factor. 


Name  of  Discusser:  De  Luca  L.,  Politecnico  Milano,  Italy 

Question: 

Impressive  3-D  calculations  were  shown  in  the  previous  presentation.  Your  technique  is  a  interesting 
mathematical  idea. 

What  specific  advantages  do  you  expect  from  the  already  available  routines: 

Accuracy  of  the  results?  Stability  of  procedure  executing  time? 

Answer: 

1  Process  of  boundary  condition  implementation  is  simplified 

2  Independent  of  used  mesh,  especially  refinement  and  type 

3  Gain  in  interpolation  speed 

4  Easy  adjustments  for  changes  in  geometry 

5  Reproduction  of  results  is  easy  because  all  important  information  is  stored  in  one  place,  namely  the 
geometric  model. 


Name  of  Discusser:  B.  Simon,  MTU  Aero  Engines  Munich 

Question: 

You  have  stored  your  boundary  condition  on  3-D  surfaces.  Does  that  mean  that  boundary  conditions  are  re¬ 
stricted  on  3  types? 

Answer? 

No,  each  surface  in  the  geometry  has  a  list  of  virtual  boundary  conditions  linked.  This  list  can  be  as  long  as 
needed  and  each  virtual  boundary  condition  contains  one  type  of  information  (e.g.  pressure,  temperature, 
heat  transfer  coeff.) 


